Universal correlations of one-dimensional interacting electrons in the gas phase 
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We consider dynamical correlation functions of short range interacting electrons in one dimension 
at finite temperature. Below a critical value of the chemical potential there is no Fermi surface 
anymore, and the system can no longer be described as a Luttinger liquid. Its low temperature 
thermodynamics is that of an ideal gas. We identify the impenetrable electron gas model as a 
universal model for the gas phase and present exact and explicit expressions for the asymptotics of 
correlation functions at small temperatures, in the presence of a magnetic field. 

PACS: 71.27.+a, 71.10.Ca, 71.10.Fd, 71.10.Pm 

There are two different phases of a one-dimensional 
system of interacting electrons. In order to explain the 
difference, let us first consider zero temperature, T = 0: 

At zero temperature there exists a critical value iic of 
the chemical potential. If /z > /Xc, the density of electrons 
in the ground state is positive. This is the familiar Lut- 
tinger liquid phase |l|. Correlations are dominated by 
fluctuations around the Fermi surface. Conformal field 
theory ^ describes their power law decay. 

If /i < fic, the density at T = is zero. This is a 
special case of the electron gas phase, which will become 
less trivial at positive temperatures. 

What happens at small positive temperatures, T > 0? 

For fi > fic the Luttinger liquid phase persists. Now 
correlations decay exponentially. In order to describe 
correlations at small temperature one has to employ a 
conformal mapping from the complex plane to a strip of 
finite width. As a result the rate of exponential decay is 
defined by conformal dimensions. 

Similarly, the gas phase persists for fi < ij,c and T > 0. 
The density of electrons is now positive, but typically ex- 
ponentially small at small temperatures. The ideal gas 
law holds. This suggested the name 'gas phase'. Corre- 
lations in the gas phase behave essentially different com- 
pared to those in the Luttinger liquid phase. This is the 
subject of this paper. 

We shall argue below that the low temperature be- 
haviour of correlation functions in the gas phase is uni- 
versal for different models. The universal model which 
determines the asymptotics of correlation functions of 
one-dimensional electrons with short-range repulsive in- 
teraction is the impenetrable electron gas, which is the 
infinite coupling limit of the electron gas with repulsive 
delta interaction. The universality class of the impene- 
trable electron gas comprises a large number of physically 
interesting systems, most notably the Hubbard model 
and its non-solvable generalizations. 

The impenetrable electron gas model is solvable by 
Bethe ansatz P,0. Based on the Bethe ansatz solution a 
determinant representation for the dynamical two-point 
Green functions was obtained in B, and, subsequently, 
their asymptotics was calculated in |^J^. The case of 



static correlations was treated earlier in 



I. THE HUBBARD MODEL IN THE GAS PHASE 

In order to get a better understanding of the gas phase 
let us start with the important example of the Hubbard 
model, 

L L 

L 

-^^(nji +nj^) . (1) 

Here the canonical Fermi operators cj^, Cj^a- are creation 
and annihilation operators of electrons at site j of a one- 
dimensional periodic lattice of length L, and nj_| and 
rij^l are the corresponding particle number operators. U 
is the strength of the Hubbard repulsion, fi is the chem- 
ical potential. 

The Hubbard model was solved Q by means of the 
nested Bethe ansatz. This allows us to test our ideas 
about the gas phase quantitatively. The energy levels for 
the N electron system are 



N 



E^ 2^(1- cos fcj)- (/i + 2)7V , 



(2) 



i=i 



where the charge momenta kj are solutions of the Lieb- 
Wu equations P] . Clearly the first term on the right hand 
side of (||) is non-negative. Hence, if /x < /ic = —2, the 
energy of all eigenstates becomes non-negative, and the 
absolute ground state is the empty lattice. For /i > —2, 
on the other hand, the energy can be lowered by fill- 
ing states with small fc's. Since kj^i — kj ~ 1/i, this 
leads to a finite density of electrons in the ground state 
as, L -^ oo. We conclude that the Hubbard model is 
in the gas phase for ^ < — 2 and in the Luttinger liq- 
uid phase else. The asymptotics of correlation functions 
of the Hubbard model in the Luttinger liquid phase was 
obtained in fp3. 



The thermodynamics of the Hubbard model was first 
considered by Takahashi [Q. He expressed the Gibbs 
free energy cu = —P {P pressure) in terms of the dressed 
energies K{k), e„(A), e^(A), of elementary excitations at 
finite temperature. K{k) is the dressed energy of parti- 
cle (or hole) excitations, £„(A) describes spin excitations 
and £^(A) fc-A strings. All k-A strings are gapped Q. 
They do not contribute to the low temperature thermo- 
dynamic properties of the Hubbard model [^ and drop 
out of the equation for the pressure, which simplifies to 



P^L 

2-K 



dk In 1 



(3) 



Similarly, the integral equations for the dressed energies 
at low temperatures become 



K{k) — — /i — 2cosfc 

oc 

-T^([n]ln( 

n=l 

InM -l-eT^ j = 



l + e" 



(sin A:) 



(4) 



dk cos k a„ (A — sin A;) In f 1 + e ^ 
^(A„,„ln(l + e-*))(A), (5) 



(6) 



where n = 1, 2, 3, ... in equation (0), and 

(\\ - nU/An 
''"^^" (n;7/4)2 + A2 ■ 

[n] and Anm, are integral operators defined by 

([0]/)(A) =. /(A) , (7) 

(H/)(A) = / dA' a„(A - A')/(A') , n = 1, 2, . . . (8) 



inin{n,7n} 

A„™= ^([|n-m| + 2(j-l)]-|-[|n-m| + 2j]). (9) 

The gas phase is characterized by the absence of a Fermi 
surface for K{k). Thus K{k) is positive in the zero temper- 
ature limit, and the first term on the right hand side of (H) 
becomes exponentially small in T. Dropping this term, 
the equations (o) decouple from (m. Since the equations 
become independent of A, it is not hard to solve them. 
The solution is the same as in the infinite coupling limit 
U ^ oo (cf. e.g. 0), exp{(r„(A)/T} = n{n + 2). Inserting 
this solution into (H) we obtain 



K{k) = -fi-2cosk-Tln2 



(10) 



Our initial assumption that lim7^_^o ^(k) > for all k is 
self-consistent, if /i -I- 2 < 0, which is precisely the condi- 
tion for being in the gas phase stated above. With (HG) 
the low temperature expression for the pressure becomes 



T I ' / M + 2C03 fc 

P= — l dk \nll + 2e — i 

27r 



T n + 2 

e^ , (11) 

TT 



and we see that density D — dP/dfi and pressure P are 
related by the ideal gas law, 



P^TD . 



(12) 



There are two important lessons to learn from our sim- 
ple calculation. First, the low temperature limit in the 
gas phase works the same way as the strong coupling 
limit at finite temperatures. Second, the low tempera- 
ture Gibbs free energy uj = — P in the gas phase shows 
no signature of the lattice. It is the same as for the im- 
penetrable electron gas (see below) , which is a continuum 
model. This fits well with our intuitive understanding of 
the gas phase: The mean free path of the electrons is 
large compared to the lattice spacing A, their momen- 
tum is small. Their kinetic energy is of the order T. 
Hence, the effective repulsion is large for T -^ U. 



II. SCALING 

The above arguments suggest that an effective Hamil- 
tonian that describes the Hubbard model in the gas phase 
can be obtained as the continuum limit of the Hubbard 
Hamiltonian. Let us introduce the continuous length 
£ = AL and coordinates x — An connected with the 
nth lattice site. In the continuum limit, A — > for fixed 
£, we obtain canonical field operators \I'o-(x) for electrons 
of spin a as 



Let us perform the rescaling 

Th = A''T , mh + 2 = AV 
tH = t/A^ , Bh - A^B , 



(13) 



kn ^ Ak , 



(14) 



where k denotes the momentum, t the time and B the 
magnetic field, which we shall incorporate below. The 
index '7?' refers to the Hubbard model. Then, in the 
limit A — + 0, we find 



Hh/Th = H/T . 



(15) 



Here H is the Hamiltonian for continuous electrons with 
delta interaction, 



H 



e/2 



dx< 



{(a,*+(:E))9,*„(a;) 
^'i+{x)^+{x)^l{x)^^{x)~ ^i^+{x)^o.{x)}. (16) 



U 



Note that the coupling ci = U/A of the continuum 
model goes to infinity! This is a peculiarity of the one- 
dimensional system. The effective interaction in the low 



density phase becomes large. Similar scaling arguments 
lead to an effective coupling C2 = U in two dimensions 
and to C3 — AC/ in three dimensions, i.e. unlike one- 
dimensional electrons three-dimensional electrons in the 
gas phase are free. 



is indeed the universal model (for small T) for the gas 
phase of one-dimensional lattice electrons with repulsive 
short-range interaction. 



IV. IMPENETRABLE ELECTRONS 



III. UNIVERSALITY 

What happens to more general Hamiltonians in the 
continuum limit? Let us consider Hamiltonians of the 
form Hq — Hh + V, where Hh is the Hubbard Hamilto- 
nian and V contains additional short range interactions. 
We shall assume that ^ is a sum of local terms Vj which 
preserve the particle number. Then Vj contains as many 
creation as annihilation operators, and the number of 
field operators in Vj is even. We shall further assume 
that Vj is hermitian and space parity invariant. 



According to equation (^^ every field Cj^cr on the 
lattice contributes a factor of A^/^ in the continuum 
limit. One factor of A is absorbed by the volume el- 
ement dx = A, when turning from summation to in- 
tegration. Thus, if Vj contains 8 or more fields, then 
V ^ A^ and V/Th vanishes. If Vj contains 6 fields, 
then at least two of the creation operators and two of 
the annihilation operators must belong to different lat- 
tice sites, since otherwise Vj = 0. A typical term is e.g. 
Vj — c~''|C^,c^j^ |Cj+i.|Cj,|Cj_|. In the continuum limit 

we have Cj+i,t = A^/^^^{x) + A^/'^d^'$^{x) + 0{A^/^). 
Hence, the leading term vanishes due to the Pauli princi- 
ple. The next to leading term acquires an additional pow- 
er of A. We conclude that V ~ A^ and thus V/Th -^ 0. 
If Vj contains 4 fields, then 



V ~ AH+{x)'^+{x)'^i{x)^^{x) + 0{A^) 



(17) 



Here the first term on the right hand side is the density- 
density interaction of the electron gas. In order to arrive 
at the impenetrable electron gas model the coefficient in 
front of this term has to be positive. Note that there 
are no terms of the order of A^ on the right hand side 
of (O) and thus no other terms than the first one in the 
continuum limit. Terms of the order of A'^ would contain 
precisely one spatial derivative. They are ruled out, since 
they would break space parity. 

Considering the case, when Vj contains 2 fields, we 
find, except for the kinetic energy and the chemical po- 
tential term, terms which correspond to coupling to an 
external magnetic field Bh- For these terms to be finite 
in the continuum limit we have to rescale the magnetic 
field as Bh = A^B (cf. equation (0)). 

Our considerations show that the impenetrable elec- 
tron gas model with magnetic field. 



Hb=H + B 



ill 



dx *+(a:X^vI,^(x) , (18) 



The impenetrable electron gas is the infinite coupling 
limit of the electron gas with repulsive delta interaction 
(A ^ 00 in (Iq)), which was the first model solved by 
nested Bethe ansatz p] . The pressure of the system as a 
function of T, /i and B is known explicitly 0], 



P= — I dfc In I 1 + e 



(19) 



and may serve as thermodynamic potential. The expres- 
sion ( |l9| ) is formally the same as for a gas of free spinless 
fermions with effective (temperature dependent) chem- 
ical potential fieff = M + nn(2cosh_B/r). Hence the 
Fermi surface vanishes for limT-»o /^e// = /^ + 1^1 < 0. 
The finite temperature correlation functions of the im- 
penetrable electron gas depend crucially on the sign of 
fj-eff- This allows us to define the gas phase at finite 
temperature by the condition fJ-eff < 0, which is also 
sufficient for deriving the ideal gas law (|2|) from the low 
temperature limit of (p^). Note that for zero magnetic 
field and small temperature equation (n9) coincides with 
therhsof (|ll|). 

The time and temperature dependent (two-point) 
Green functions are defined as 



G||(x, tj 
G||(a;, t) 



tr(e-«/^*t(x,t)*+(0,0) 



tr (e-H/T) 
tr(e-^/^*+(x,t)*|(0,0) 



tr (e~"/n 



(20) 
(21) 



For the impenetrable electron gas these correlation func- 
tions were represented as determinants of Fredholm in- 
tegral operators in |B|. The determinant representation 
provides a powerful tool to study their properties analyt- 
ically. 

The short distance asymptotics can be obtained 'per- 
turbatively' by expanding the determinant representa- 
tion for small x and t and fixed ratio ko — x/2t, 

G+^{x,t) = ^-L + ^{kl+fi- B)Vt 
2v7i'i ^v"" 

i-^)f-(>.^)f.o<„,,. 



Gf^ (x, t) = Dt + ^t it + O (i3/2^ 



(23) 



Here 7 = 1 -h e^/^, D-^ = dP/d{fi - B) is the density of 
up-spin electrons and 



Er — — — 
2tt 



dk- 



(fc2 - fj, + B)e~~ 
2cosh{B/T) + e^ 



(24) 



may be interpreted as the 'energy' of the up-spin elec- 
trons. In the gas phase at low temperatures the density D 
becomes exponentially small and the terms proportional 
to D may be neglected in (|2^). Then, to lowest order in 
D, (|23) and ( P3| ) are the same as for free fermions. 

In [|6|J7| the determinant representation was used to 
derive a nonlinear classical differential equation, which 
drives the correlation functions. This equation is closely 
related to the quantum Hamiltonian (p^). It is the sep- 
arated nonlinear Schrodinger equation. Together with 
a corresponding Riemann-Hilbert problem it determines 
the large-time, long-distance asymptotics of the correla- 
tors (pO|), (|2l|). In 1^,0 the asymptotics x,t ^ oo was 
calculated for fixed ratio ko = x/2t. The crucial parame- 
ter for the asymptotics is the average number of particles 
xD in the interval [0, x]. li x is large but xD ^ 1 (i.e. T 
small), an electron propagates freely from to x, and the 
correlation functions (^) , (|2^) are those of free fermions. 



G+{x,t) 
G-{x,t) 



2V^ 



2V^ 



f-ipitCM-s) 



r 



^e-i*(M-s)g- 



(25) 
(26) 



The true asymptotic region is characterized by a large 
number xD of particles in the interval [0, a;], specifically, 
xD > z-i, where z, = (r3/4e-'=o/2T)/(2^i/4fc3/2). if 
the latter condition is satisfied, the correlation functions 
decay due to multiple scattering. The cases B > Q and 
B <Q have to be treated separately. For B > Q there is 
a critical line x = At^fB in the x-t plane, which separates 
it into a time and a space like regime. The asymptotics 
(for small T) in these respective regimes is: 

Time like regime {x < 4t\/B): 



G+^{x,t)=G+{x,t) 



r 



il/(Ze), 



E_D| 



^JA■nZcxD<^ 



G^-^{x, t) 



G-Jx,t) 



fil^iZa}^ 



-xD, 



y/AnZcxDi 



where 



iy{zc) = - 



2D^kl'''e-^"/^^ 

^l/4y5/4 



D, = \J^e'^^+B)/T 



Di^ = dP/d{fi + B) is the low temperature expression for 
the density of down-spin electrons. 

Space like regime (x > 4t^/B): 



G+(x,t)=G+(x,t)t-''^(^ ')e-^-^i , 
ijT -1- -^. yx ^ t J =^ ^ f yX J tj t 



where 



Hi-') 



i'^h )p-xDi 



„(3B+M-feo)/^ 



2n 



(31) 

(32) 



(33) 



For B < there is no distinction between time and space 
like regimes. The asymptotics is given by (|3l]), (p^). 

Equations (|2^), ( pq ) and (|3l|), ( p2| ) are asymptotic ex- 
pansions in t consisting of an exponential factor, a power 
law factor and a constant factor. Note that the method 
employed in [Q allows for a systematic calculation of the 
next, subleading orders. 

The leading exponential factor in (|27|), ( [2q ) and (|3l|), 
( P2| ) has a clear physical interpretation: Because of the 
specific form of the infinite repulsion in (|l6|) , up-spin elec- 
trons are only scattered by down-spin electrons. This is 
reflected in the fact that the correlation length is l/D^. 
The exponential decay means that, due to the strong 
interaction, an up-spin electron is confined by the gas 
of surrounding down-spin electrons. Thus we are facing 
an interesting situation: Although at small distances the 
electrons look like free fermions, they are confined on a 
macroscopic scale set by the density D^ . 
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